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^! Abstract 
< 

, ' We consider a continuum of phase oscillators on the circle interacting through an impulsive 

'NT ; 

instantaneous coupling. In contrast with previous studies on related pulse-coupled models, the 

^ . stability results obtained in the continuum limit are global. For the nonlinear transport equation 

' governing the evolution of the oscillators, we propose (under technical assumptions) a global 

' Lyapunov function which is induced by a total variation distance between quantile densities. 

I The monotone time evolution of the Lyapunov function completely characterizes the dichotomic 

behavior of the oscillators: either the oscillators converge in finite time to a synchronous state 

or they asymptotically converge to an asynchronous state uniformly spread on the circle. The 

results of the present paper apply to popular phase oscillators models (e.g. the well-known leaky 

integrate-and-fire model) and draw a strong parallel between the analysis of finite and infinite 

Tlj" , populations. In addition, they provide a novel approach for the (global) analysis of pulse-coupled 

(N . oscillators. 
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I. Introduction 

Networks of interacting agents are omnipresent in natural [B], [21], [3Z1 as well as in 
artificial systems [12], [21], [SI]- In spite of their apparent simplicity, they may exhibit 
rich and complex ensemble behaviors [32] and have led to intense research during the last 
few decades. In this context, coupled phase oscillators are generic models of paramount 
importance when studying the collective behaviors of a large collection of systems |38j . 

Phase oscillators appear as reductions or approximations of (realistic) dynamical oscillator 
models. They are obtained through the computation of a phase response curve (PRC) |3Uj . 
|38] . which characterizes the phase sensitivity of an oscillator to an external perturbation, 
such as the influence of the neighboring oscillators in the network. Since phase oscillators 
are characterized by a one- dimensional state-space 5*^(0, 27r), they are more amenable to a 
formal mathematical study of the collective behaviors even though the nonlinear interactions 
between oscillators often yield mathematical puzzles [I6], [36] . 

Within the network, oscillators interact through a nonlinear coupling. In most studied 
models, the coupling has a permanent influence on the network. However, in many situations 
encountered in biology or physics, the oscillators influence the network only during a tiny 
fraction of their cycle (e.g. yeast cell dynamics [2]). It is particularly so when the interconnec- 
tion between the agents consists in the emission of fast pulses (spiking neurons [TU] , cardiac 
pacemaker cells [26], earthquakes dynamics [23], etc.). In this paper, we consider the limit of 
an impulsive and instantaneous coupling, where a (pulse-coupled) oscillator interacts with 
the network only when its phase is equal to a given value. When considering the popular 
leaky integrate-and-fire oscillators, this model corresponds to Peskin's model [26], [T7] . 

For a finite number of identical oscillators, previous studies show that the global behavior 
of pulse-coupled oscillators is dichotomic ([IS], [20]): the oscillators converge either to a 
synchronized state or to an ant i- synchronized state. In the present paper, we extend the 
result to infinite populations, showing that the global behavior of infinite populations is 
also dichotomic and thereby highlighting the perfect parallel between finite and infinite 
populations. 

Several earlier studies have provided local stability results for infinite populations of pulse- 
coupled oscillators (see e.g. [T], [5], [?], [T7j, [33]). In contrast, we present in this paper global 
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stability results for infinite populations of monotone oscillators (including leaky integrate- 
and-fire oscillators). To this end, we introduce a Lyapunov function which is induced by a 
norm and which has the interpretation of a total variation distance between (quantile) 
density functions. Modulo technical conditions detailed in the paper, we show that the time 
evolution of the proposed Lyapunov function is governed by the derivative of the PRC, a 
result that leads to a global convergence analysis for monotone PRC's. 

Beyond the analysis of monotone pulse-coupled oscillators, the theory developed in the 
paper leads to general results on nonlinear partial differential equations (PDE). The analysis 
focuses on transport equations with a monotone dynamics (derived from a monotone PRC) 
and provides existence, uniqueness, and global stability results for the stationary solution of 
the PDE. In particular, the use of a total variation distance as a strict Lyapunov function 
seems novel and specific to the impulsive nature of the coupling (the results in [22], [33] 
suggest that total variation distance is of little use for systems of conservation laws). In this 
context, the result could potentially open new avenues to connect the monotonicity property 
of an infinite dimensional system to its stability properties. 

The paper is organized as follows. In Section [Tll we derive the transport PDE for the 
continuum of pulse-coupled oscillators from the original model of pulse-coupled integrate- 
and-fire oscillators (Peskin's model). Section IIIII presents numerical experiments showing 
that the continuum model is characterized by a dichotomic behavior. In addition, existence 
and uniqueness results are obtained for the stationary solution of the PDE. In Section IIVI a 
strict Lyapunov function is proposed, which is inspired from our previous work [12] on finite 
populations. In Section |Vl we perform the convergence analysis of populations of monotone 
oscillators. The parallel between finite and infinite populations, as well as some extensions 
of the model, are discussed in Section IVI[ Finally, the paper closes with some concluding 
remarks in Section IVIIi 

II. A PHASE DENSITY EQUATION FOR PULSE-COUPLED OSCILLATORS 

In this section, we introduce models of (monotone) integrate-and-fire oscillators with an 
(instantaneous) impulsive coupling. The oscillators are equivalent to phase oscillators and, 
in the continuum limit, evolve according to a phase density equation. The derivation of the 
corresponding nonlinear PDE is standard and similar developments are found in [T], [S], |17j . 
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A. Pulse-coupled integrate- and- fire oscillators 

We consider models of integrate-and-fire oscillators [13]. An integrate-and-fire oscillator 
is described by a scalar state variable x, which monotonically increases between the two 
thresholds x and x according to the dynamics x = F{x), F > 0. When the oscillator reaches 
the upper threshold x, it is reset to the lower threshold x (it is said to fire). 

In Peskin proposed to study the behavior of integrate-and-fire oscillators inter- 
acting through an impulsive coupling. Whenever an oscillator fires, it sends out a pulse 
which instantaneously increments the state of all other oscillators by a constant value K/N, 
where K is the coupling strength. The coupling is usually excitatory {K > 0) but may 
also be inhibitory {K < 0). The dynamics of a pulse-coupled integrate-and-fire oscillator 

G {1, . . . , A^} is then given by 

Xk = F{xk) + Uk(t) (1) 

with the coupling 

jy' N oo 

Mt) = ^T.T.^it-t?^)- (2) 

i=l 1=0 

The Dirac functions 5 model the pulses which increment the state of oscillator k at the firing 
times tp'', that is, when an oscillator i fires. 

Peskin's model was initially proposed with the popular leaky integrate-and-fire (LIF) 
oscillators, characterized by the affine vector field F{x) = 5* — 7X, 7>0 [TS]. However, the 
results of the present paper apply to monotone oscillators, that we define as integrate-and- 
fire oscillators with a monotone increasing or monotone decreasing vector field {dF/dx > or 
dF/dx < 0). The class of monotone oscillators embraces a large variety of models, including 
the popular LIF model. 

B. Phase oscillators 

Integrate-and-fire oscillators are equivalently modeled as phase oscillators if the state 
dynamics ([T]) is turned into a phase dynamics. The phase 6 G 5*^(0,27?) is determined from 
the state x G [x,x] by rescaling in such a way that ^ = corresponds to the low threshold 
X = X — the oscillator fires at phase ^ = — and in such a way that a single (uncoupled) 
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oscillator has a constant phase velocity 6 = u, where u is the natural frequency of the 
oscillator. This leads to the state-phase relation 

Under the influence of the coupling Ufc(t), the state dynamics ([T]) corresponds to the phase 
dynamics 

ek = u + z{ek)uk{t)^viek,t), (4) 

where the function Z G C^([0,27r]) is the infinitesimal phase response curve (PRC) of the 
oscillator, that is, the phase sensitivity of the oscillator to an infinitesimal perturbation [T3j, 
|3U] . For integrate-and-fire oscillators, the PRC has the closed-form expression (see 

It follows from that monotone oscillators are characterized by a monotone PRC: dF/ dx > 
Vx e [x, x] leads to < G [0, 27r] and dF/dx < Vx G [x, x] leads to > G [0, 27i] , 
where Z' denotes the derivative of Z with respect to the phase 6. 

C. Phase density equation 

In the limit of a large number of -> oo oscillators, the inGnite population is a continuum 
characterized by a (nonnegative, continuous) phase density function 

p{e,t) G C°([0,27r] X R+;M+) 

that satisfies the normalization 

["^^ p{d,t)dd = l Wt. 
Jo 

The quantity p{6,t)d6 is the fraction of oscillators with a phase between 6 and 6 + d6 a,t 
time t. The time evolution of the density obeys the well-known continuity equation 

§-^Pio,t) = -^[v{e,t)p{e,t)] , (6) 

where the function 

vie,t)pie,t)^jie,t) (7) 

is the (nonnegative, continuous) flux J{6,t) G C°([0,27r] x R+;M+). The quantity J{9,t)dt 
represents the fraction of oscillators flowing through phase 9 between time t and t + dt. 
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Since the phase 6 is defined on S'^(0,27r) = M mod 27r, the fiux must satisfy the boundary 
conditions 

J(0,t) = J(27r,t) = JoW Vt. (8) 

For the sake of simphcity, we use in the sequel the notation Jq to denote the boundary 
fiux ([H]). Since the oscillators fire at phase 6* = 0, Jo(^) is also called the &nng rate of the 
oscillators. 

D. Continuous impulsive coupling 

The impulsive coupling originally defined for finite populations is extended to infinite 
populations as follows. Since the coupling strength K/N is inversely proportional to the 
number of oscillators, the coupling does not increase as the number of oscillators grows: the 
constant K corresponds to the net infiuence of the whole population through the coupling, 
regardless of the number of oscillators. In the limit — ?■ oo, the firing of each oscillator of 
the infinite population produces an infinitesimal spike of size K/N — )■ 0. 

The impulsive coupling is best expressed in terms of the fiux Jq. For a finite population, the 
oscillators crossing 6* = at times induce a (discontinuous) fiux Jo(t) = J2j J2i ^{t — ^F^)- 
In the limit N oo, the infiuence of a single oscillator is negligible, so that comparing the 
fiux Jo with the coupling ([2]) yields the coupling 

u{t)=KJo{t). (9) 

Roughly, for infinite populations, the impulsive coupling is identical for all the oscillators 
and proportional to the firing rate Jq. In addition, the coupling is a continuous-time function 
interpreted as an infinite sum of infinitesimal spikes that result from the uninterrupted firings 
of the continuum. 

With the continuous impulsive coupling the phase dynamics is rewritten as 

e = v{e,t) = uj + z{e)KJo{t). (lo) 

The final PDE for the continuum of pulse-coupled oscillators is derived as follows. At 6* = 0, 
([7]) and (ITU]1 yield the relationship 

Jo{t) = v{0,t)p{0,t) = [ij + KZ{0)Jo{t)]p{0,t) (11) 
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and the flux at the boundary is exphcitly given by 

ujp{0,t) 



l-KZ{0)p{0,t) ■ 

The continuity equation (jH]) thus leads to the nonlinear PDE for the density 
dp{e,t) dp{e,t) KLop{0,t) d 



-Ld- 



dt 36 I - K Z{Q) p{Q,t)de 

The boundary condition ([8]) is expressed in terms of density as 

p(0,t) p(27r,t) / Jo(t) 



[z{e)p{e,t)] . (12) 



(13) 



1- K Z{0)p{0,t) 1- K Z{2tt) p{2n,t) \ u J ' 
The reader will notice that, in contrast to the periodicity condition on the flux, no periodicity 
is assumed on the density p{0,t). [In particular, p(0,t) p{2iT,t) if Z{0) ^ Z{2'k).] 

The density of the pulse-coupled oscillators evolves according to the nonlinear PDE 0121) 
with the boundary condition (fT3l) . The PDE is studied in detail in the rest of the paper, 
with a particular attention to the case of a monotone PRC. Observe that a monotone PRC 
implies that Z(0) ^ Z(27r). 

III. A DICHOTOMIC BEHAVIOR 

Finite populations of pulse-coupled monotone oscillators exhibit a dichotomic behavior: 
they converge either toward a synchronized state or toward an anti-synchronized state 
(see [in], [20] )• Similarly, inflnite populations are characterized by a dichotomic asymptotic 
behavior, that depends on the coupling sign {K > or A' < 0) and on the derivative Z' (or 
equivalently dF/dx). This remarkable behavior is described in the present section through 
numerical experiments and intuitive arguments, motivating the theoretical global analysis 
in the next sections. The actual proof of the dichotomic behavior is postponed to Section 

m 



A. Asymptotic behavior 

Without coupling [K = 0), the last term of f[T^ disappears and the PDE is a standard 
transport equation. Its solution is a rigid translation of the initial density p{9,0) = po{0) 
with a constant velocity u, that is, a traveling wave p{0,t) = po{{0 — ut) mod 2ti). In this 
case, any solution is periodic (with period 27r/a;) and the system is marginally stable. 
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When the oscillators are coupled, the last term of fll2l) modifies the transport equation. 
Under the infiuence of the coupling, the velocity depends on both time and phase and the 
density is thereby "stretched" or "compressed". This is illustrated when computing the total 
time derivative along a characteristic curve A{t) defined by A = v{A{t),t), that is 

^ + I^H^W'O = f = Mt)KZ'{A{t)) , (14) 

where f lT^ and (IT^ have been used. The total derivative shows that the density is modified 
on a characteristic curve whenever the PRC Z is not constant. In addition to the rigid 
translation, the density undergoes a nonlinear transformation, possibly leading to asymp- 
totic convergence to a particular density function corresponding to a particular stationary 
organization of the oscillators. 

The total derivative (IT^ gives clear insight that the sign of the derivative K Z' is of 
primary importance. In fact, the sign of K Z' will enforce a dichotomic behavior. The 
condition K Z'{9) < (or KdF/dx > 0) will be shown to enforce convergence to a uniform 
flux J{6,t) = J* on 5^(0, 27r). This situation, corresponding to the maximal spreading of 
the oscillators on the circle, is called the asynchronous state [T], [53] (Figured]). In contrast, 
the reverse condition K Z'{9) > (or KdF/dx < 0) will be shown to enforce convergence 
to a delta-like flux (Figure [21). This situation, characterized by the synchronization of all the 
oscillators, is the synchronous state. 




% 9/2Tr 1 5 t 10 15 

(a) (b) 

Figure 1. When K Z' < (or K dF/dx > 0), the solution converges to the asynchronous state. With the monotone 
LIF dynamics x — 2.1 — 2x, x £ [0, 1] and with an inhibitory coupHng K = —0.1 < 0, the function K Z is monotone 
decreasing, (a) The density converges to a stationary solution p* and (b) the flux Jo{t) tends to a constant value 
J* 0.53. 
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Figure 2. When K Z' > (or K dF/dx < 0), the solution converges to the synchronous state. With the same LIF 
dynamics as in Fig. [1] but with an excitatory couphng K = 0.1 > 0, the function K Z is monotone increasing, (a) 
The density converges to a synchronous solution and (b) the flux Jo{t) tends to a Dirac function. 



We remark that reversing the couphng sign {K > or < 0) has the same effect as 
reversing the monotonicity of Z (Z' < or > 0). 

The dichotomic asymptotic behavior is in agreement with the equivalent dichotomic 
behavior of finite populations: the synchronous and asynchronous states are the exact analog 
of the synchronous and splay states observed for finite populations. This is obvious in the case 
of synchronization, that occurs both for finite and infinite populations when K dF/dx < 
|2U] . Asynchronous state and splay state are also equivalent behaviors, that both occur when 
K dF/dx > ^19\. They are both anti-synchronized behaviors, for which the oscillators 
evenly spread over the circle 5^(0,27?). 

B. Stationary asynchronous state 

We will now characterize the stationary solution of the PDE (j6]), i.e. the asynchronous 
state corresponding to a constant flux J{0,t) = J*. From ([7]), the stationary density must 
satisfy 

''^"'- .Jmj' - '''' 

The stationary asynchronous state thus exists if there exists a value J* > so that the 
stationary density is nonnegative and bounded 
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and normalized 

The condition (I16p and the continuity of Z imply that the stationary solution is continuous. 

It is noticeable that for finite populations, the equivalent stationary splay state is a phase- 
locked configuration: at each firing, the N oscillators are characterized by constant phases 
0* = {61, ■ ■ ■ ,6^_^,d^ = 2it) and fire at a constant rate. In particular, the value ^tv-i is 
related to the firing rate and is well-approximated by the stationary flux J* when ^ 1. 



The following proposition gives necessary and sufficient conditions to ensure the existence 
and uniqueness of a stationary solution p* G C°([0, 27r]; M"*"). 

Proposition 1. A stationary flux J* > satisfying the conditions l[T5\) and l[T7\) exists if 
and only if the inequality 

lim / -rrTTTTT, dO > 1 (18) 

^^^;7o KZ{e) + s ^ ' 

is satisfied with 

ifKZ{e)>0 V^e[0,27r], 
miiiig^[o^2n] (^K Z (9)^ otherwise. 
Moreover, the solution is unique when it exists. 

Proof: Inequality flTB]) implies that the velocity uj + K Z{6)J* is strictly positive, so that 

s>r 




The function 

W{.J) = / TTTTTT^ de 



/ ^d9>0 WJeJ. 

Jo \lo + KZ(B)JV 



uj + KZ{e)J 

satisfies W{0) = 0, is continuous on J', and is strictly increasing on J' since 

dW _ f^^ CO 
U ~ Jo [uj + KZ{9)JY 

As a consequence, the equation W{J) = 1, which is equivalent to condition fll7l) . has a 

(unique) solution J* G if and only if 

lim W{co/s) = lim / — — d9 > 1 , 

^ ' ' ^^^Jo KZ(9) + s 



s>r s>r 
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which concludes the proof. ■ 
Proposition [Uimphes that the couphng constant K must be bounded. For integrate-and- 
fire oscillators, the bounds on the coupling constant are computed analytically and are given 
in the following corollary. 

Corollary 1. For integrate-and-fire dynamics x = F{x), a stationary flux J* > fulfilling 
the conditions ( fl^) and fTl\} exists if the coupling constant satisfies 

rx g 

lim / ——dx<K<x — x, (19) 



S<-Fmin 

x) 



with Fmin = mina;e[^,3;] {f{ 

Proof: If the coupling is excitatory {K > 0), it follows from that K Z{6) > 
V6' G [0,27r], so that r = 0. Then, the condition (fT8l) of Proposition [1] can be rewritten as 

- 1 ,,^/-fM^,,= ri,,>i 



/o KZ{e) Jo Ku Jx K 

given ([3]), or equivalently 

K <x-x. (20) 

If the coupling is inhibitory [K < 0), it follows from that r = —Ku/F^^jy. Then, the 
condition ( fTSi) of Proposition [1] is rewritten as 

lim / -— dO = lim / — ^^dO = lim / uj-^dx > 1 

and we obtain the lower bound (HM on K. ■ 
Condition (1201) is easy to interpret in the case of finite populations: it is a necessary 
condition for the existence of a stationary phase-locked configuration. Since the average 
state difference between oscillators is (x — x)/N, a coupling strength K >x — x yields an 
"avalanche" phenomenon igniting a chain reaction of firings. In this situation, a phase-locked 
behavior of distinct oscillators cannot exist. 

IV. A STRICT Lyapunov function induced by the total variation distance 

Lyapunov analysis is a classical approach to study the stability of nonlinear PDE's (see 
e.g. [7]). In this section, we extend our previous results obtained for finite populations [TH] 
to construct a strict Lyapunov function for the PDE ( fT2l) . The Lyapunov function, inspired 
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by the 1-norm introduced in [T^], is a distance interpreted as the total variation distance 
between quantile densities. 

A. Quantile density 

The description of the infinite population through the density p{6) is not equivalent to 
the description of a finite population through a vector = (6*1, ■■ ■ , Oj^i). While the former 
corresponds to the "amount" of oscillators as a function of the phase, the latter corresponds 
to the phase as a function of the oscillator index. To establish an equivalence between finite 
and infinite populations, we introduce an index for infinite populations of oscillators and 
use the concept of quantile function. 

For infinite populations, the oscillators can be continuously labeled on the interval [0, 1] 
and an oscillator index if G [0, 1] is defined as follows. Given a density function p : [0, 27r] i— > M"*", 
the cumulative density function P{9) : [0,27r] i— )• [0, 1], defined as 

P{9)= f%{s)ds, 
Jo 

attributes an index ip = P{9) G [0, 1] to each oscillator with phase 9. In particular, an index 
(/:> = (resp. (p = 1) is attributed to the oscillator at phase 6' = (resp. 9 = 27r). 

Next, to complete the equivalent description of infinite populations, we introduce the 
quantile function (widely used in statistics [25]): the quantile function Q : [0, 1] t— ^ [0,27r] is 
the inverse cumulative density function, that is, 

Q{if) = p-\^) = inmP{9) > ^} . 

The (continuous) quantile function is equivalent to the (discrete) description of finite 
populations. For a finite number of distinct oscillators, at each firing of an oscillator 
{9n = 27r), the remaining phases 9k are the A^-quantiles 9^ = Qk^\ that is 

= (^1, ■ ■ ■ , 27r) = {Qf\ ■ ■ ■ , gS?2i, 27r) . (21) 

When the number of oscillators tends to infinity, the A^-quantiles are replaced by the 
continuous quantile function Q. (Roughly speaking, the quantile function plays the role 
of the vector with an infinity of components.) 
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As the analog of the density p, the quantile density function [23], also called sparsity 
function, is the function q : [0, 1] i— )• that satisfies (see Figure [3]) 

Qm = / q{s)ds. 
Jo 



The quantile density function, which is the derivative of the quantile function, expresses the 
increase of phase per unit increase of oscillator index. The density function is linked to the 
quantile density function by the relationship 




Figure 3. The density function p{9) (left) has a cumulative density P{9) (center). The quantile function 
Q(<p) = P^^if) is the cumulative density function of the quantile density function q{'p) (right). 

The reader will notice that, as the oscillators density p{6, t) depends on time in the model 
the associated quantile function and quantile density function also depend on time and 
are then rigorously defined as the two-variable functions Q{ip,t) and q{ip,t). In addition, we 
denote the quantile function and the quantile density associated to the stationary solution 
f lTSll by Q*{(p) and q*{(p) respectively. 

B. Total variation distance 

A 1-norm introduced in our previous study [TH] on finite populations of monotone pulse- 
coupled oscillators leads to a distance (between quantile densities) in the case of infinite 
populations, a distance which can be interpreted as a total variation distance. 
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For a finite population of distinct oscillators, it is sufficient to consider only the firing 
instants, for which 6'jv = 27r. Then, the corresponding discrete system is characterized by a 
simple Lyapunov function: the distance (induced by a 1-norm) between a configuration 
and the stationary phase- locked configuration 0*. It is expressed as 

Af-2 



v(^) = - ei\ + \iOk - ^.+i) - {ei - ei_,. 



k=l 



N-l 



9 



N-l 



The phases 9k can be replaced by the iV-quantiles, according to (ETj) . and one obtains 



*{N) 



N-2 
k=l 



(TV) 



Vfe+1 



^*(TV) 
Vfe+1 



^(TV) _^*{N) 
Vat-I VTV-1 



(23) 



In the limit — t- oo, the continuous equivalent of 
between the quantile density functions: 



corresponds to the distance 



V(P) 





dQ 


dQ* 


Jo 


dip 


dip 



dp 



\q-q*\\Li Vp e C"([0,27r] X R+;M+). (24) 
One verifies that V(p) = -v^ g = g* a.e. -v^ 



The second equality is obtained through 
p = p* a.e. 

Our previous study [19] shows that, under mild conditions, quantity (l23ll decreases at the 
successive firings of the oscillators, enforcing a contraction property for the 1-norm. We claim 
that, for infinite populations, the continuous equivalent f l24p also decreases at the successive 
firings of the continuum, that is, ( IMl) decreases continuously with time. The main result of 
this paper will thus establish fl24p as a good Lyapunov function for the PDE f[T^ . 

The Lyapunov function f l24p is interpreted as a total variation distance. Indeed, the total 
variation distance between two random variables corresponds to the distance between 
the corresponding density functions (see [5] for further details). In the present case, the total 
variation distance is the sum of the maximum differences between the two quantile functions 
Q and Q*, minus the sum of the minimum differences (Figure HI). 



C. Time evolution of the Lyapunov function 

Our main result shows that the candidate Lyapunov function 
evolution provided that the PRC is monotone. 



has a monotone time 
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Figure 4. The Lyapunov function (I24p is the total variation distance between two quantile density functions. In the 
example of the figure, the distance is equal to V = 2(Aa + Ac + Ad) — 2Ab, with Aa,b.c,d > 0. 



Theorem 1. Let p{e,t) G C°([0, 27r] X W^; Mq"), with ||p||li — l.> be a strictly positive solution 
of fT^) - fTS\) . If the stationary density ( 173]) exists and if either Z"[9) > 76* G [0,27r] or 
Z"(9) < V6' G [0, 27r], then the Lyapunov function satisfies 



J(0,t)^min^ ^KZ'{9))V{p) < V(p) < ^(0,t) max^ [K Z'{9)jV{p) . (25) 
Proof: The time derivative of can be written as 

^^P^ = Jo [d^-^)dt[d^-^)'^^ = Jo i j ' ^''^ 

with G{(f,t) = dQ/dif — dQ* /dip and with the signum function sign(x) = 

Differentiating the expression 9 = Q{P{9,t),t) with respect to time t leads to 

= j^[Q{P{e,t),t)] = ^{P{9,t),t) + ^ino,t\tf-^{9,t) 



or 



ot dt dip dt p[Q{ip,t),t) 



given ( l22l) . Furthermore, 



dP do r^dJ 

-i9,t)= -£(s,t)ds = - —(s,t)ds = Mt)-Ji9,t) 



dt' ' ' Jo dt' ' ' Jo d9 

and the expression (1261) becomes (from this point on, we drop the time variable) 
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Assume that the function G{(p) has a finite number of zero crossings, that is, a finite 
number of values = < v?^^^ < ■ ■ ■ < (y?^^'^^ < (p[^''~^^^ = 1 satisfying Gi^ipf^) = with 
either a right or a left nonzero derivative. The assumption on Nc will be relaxed at the end 
of the proof. Without loss of generality, assume that G is nonnegative on even intervals and 
nonpositive on odd intervals, that is, 

'dQ dQ* 



-1 



dip dip 

Next, computing the integral in (1271) . we obtain 



> G [(^ 



Nc 



V(p) = 2E( 



(28) 



k=i p[Q{.V' 

Since it follows from (l22l) and the equality G{ip[''^) = that 



(29) 



one obtains, given ([7j), ffTUl) . and (IT^ . 



V(P) 



2E(-i) 

fc=l 

2E(-i) 

fc=l 



Jo 



P*{Q*{ip^^ 



+ i^z(g(^f)))jo - '^-Kz{Q*id''))Jo 



(30) 



The boundary condition (1131) yields a monotone relationship between the values p(0) and 
p(27r), that is, p(0) > p*(0) if and only if p(27r) > p*(27r). Apart from the case p(0) = p*(0), 
p(27r) = p*(27r), the number N^, of values ip^'^^ satisfying fl2^ . is even, owing to the continuity 
of p and Q (p* and Q*). Consequently, the terms {—l)^{u — uJq/J*) in fl5Ul) cancel each 
other. In the particular case p(0) = p*(0), is not necessarily even but it follows from (1131) 
that Jo = J*{0) = J* and the above-mentioned terms are equal to zero. As a consequence, 
one obtains 



Nc 



V(p) 



2JoE(-l)'^ 



k=l 

Nc 



{ky 



{ky 

c , 



k=l 

Nc 

= 2JoE^^'^ 



(31) 



k=l 
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where the second equahty is obtained through the mean value theorem, with e [Qi'^i'''^), Q*i^c'^'^)] 
or 6G[Q*(^f)),g(^W)]. 

It remains to consider separately each term T^'^^ in the sum ( 1311) . For the sake of simplicity, 
we first consider the case K Z" > 0. Denoting {—l)^[Q{ip^J^'>) — Q*{(p[^^)] by A^^^Q, we also 
distinguish two cases: A^'^^Q > and A^'^^Q < 0. 

Case A('^)Q > 0. One has 



rp{k) ^ KZ'{ik)A^^^Q < m^^^[KZ'{e) ]A^^'^Q (32) 



and 



= KZ'(a)AWQ > min ^(^KZ' {6) jA^^^Q. (33) 

Case A'-'^^^Q < 0. We need to consider the addition of the term T'-'^^ with the term T^^~^^ 
or T^^+^\ By ([21, one gets 

- A(^-i)Q < A^'^^Q < 0, (34) 

_A('=+i)q < A^^^Q < 0. (35) 

It follows that 

y(fc-i) ^ rj^ik) ^ K Z'{ik-i)A^^-^^Q + KZ'{ik)A'^^^Q < KZ\^k) (A^'^'^^Q + A^^^Q) . (36) 



The assumption K Z" > implies K Z\^k-i) < K Z'{$,k), with ^k-i < C,k and the above 
inequality then follows from (IMl) . In addition, also implies that the right hand in 
inequality (13^ is the multiplication of K Z' with the positive quantity A'^'^^^^Q + A^^^Q. 
Hence, can be rewritten as 

y (fc-i) ^ < ^^^^ ^ (^^ ^/(^)^ (^A(*^-i)Q + A^^'^q) . (37) 
Similarly, considering the addition of the terms T^'^^ and T^'^^^^ and using fl5S]) . one obtains 
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Next, the inequalities (1321) and (13 7p imply that the expression (13 ip can be rewritten as 

V{p) = 2JoY.T^'^ < 2^0 max A'Z'(^)U: A^Q 

k=i 0m-2n] \ / f^^^ 



Jo max (KZ'ie)) /' sign f ^ 

9G[0,27r] V J Jo \d(p 



dQ 


dQ*\ 


dQ 


dQ*\ 




dip J 




dip J 



= Jo max [KZ'(e) ]V(p) . (39) 

In the case A(^')g < 0, the two terms T^^') and T^^~^^ are considered together. The 
additional T^^~^^ itself corresponds to the case A*-^'^^^(5 > 0, and does not require to be 
associated in turn with another term. In addition, there is no boundary problem since the 
term T^^) satisfies A^^^Q > 0, given ^ and g(0) = Q*(0) = 0. 

Similarly, the inequalities (1330 and (I38p lead to 

V(p) = 2 Jo T^'^ > Jo ^min [k Z\e)\ V(p) . (40) 



In the case K Z" < 0, the inequalities ( 1371) and (l38l) are reversed, that is, the sum T^^ + 
T^^'^ has a lower bound and the sum T'^^'' +T^^^^^ has an upper bound. Hence, the inequalities 
dSni) and (HOD still hold. 

We have completed the proof assuming a finite number of zero crossings G{lp^^^) = 0, 
but the reader will notice that the bounds fl5I?]l and fHUl) do not depend on Ai'c. For an 
arbitrary p, we conclude the argument as follows: by continuity, there exists a sequence of 
Pi p such that every member of the corresponding sequence Gi{ip) — )■ G{ip) has a finite 
number of zero crossings. Because the bounds fl5I?]) and fHU]) hold independently of i, they 
also hold in the limit, which concludes the proof. ■ 

If there is no coupling or if the PRC is constant, (IT^ is a (marginally stable) standard 
transport equation and Theorem [1] implies that the Lyapunov function is constant along 
the solutions of (I12p . This is in agreement with the fact that the total variation distance 
is a conserved quantity for most of the systems of conservation laws [22]. But whereas the 
distance is constant with a standard transport equation, the distance is not constant under 
the infiuence of the coupling: the nonlinear coupling term in f[T21) . which depends on K Z', 
induces a variation (I25p of the Lyapunov function. A monotone decreasing function K Z 
implies a decreasing Lyapunov function along the solutions. This is discussed in detail in 
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the next section. 



Vus{p)= r\p-p*\de VpGC°([0,27r] x]R+;]R+). (41) 
Jo 

e is given by 

Vbis = / sign(p - p*)^de = y: (-1)' [Aoi'-'^) - M'^)] , (42) 

Jo at tz^ ^ J 



The above result emphasizes the importance of considering (i) quantile densities instead 
of densities and (ii) a distance (total variation distance) instead of a distance. The 
importance of these two points is illustrated in the two following paragraphs. 

Quantile density vs. density function. The Lyapunov function 0241) is induced by the 
total variation distance between quantile densities. An alternative choice would be the total 
variation distance between density functions, that is, 

10 

The time derivative is given by 

k=l 

where we used (E]) and introduced the notation analog to that is, the values 9[''^ 

{k = 1, . . . , Nc) satisfying 

p{ei'\t) = p*{e['\t). (43) 

A simple argument shows that the Lyapunov function cannot be strictly decreasing along 
the solutions of ( fT2l) . Indeed, for any density satisfying p(0) = p*(0), one must also have 
Jo = J*(0) = J* since it follows from f[T5]) that there is a bijection between the values p(0) 
and Jo- But then, and (|43]) imply that J(6'('')) = J* {6^'^) = J* for all k and the derivative 
( H2l) leads to Vbis = 0. Even though one actually shows that V < 0, La Salle principle cannot 
be used to prove the stability, since it is not obvious to prove the precompactness of the 
trajectories. 

This argument shows that a direct application of the total variation distance on density 
functions does not lead to a good candidate Lyapunov function for the PDE (fT21) . A key 
point is to apply the total variation distance on quantile functions instead. 

L-*^ distance vs. distance. The Lyapunov function (1241) is induced by a distance. An 
alternative choice would be a Lyapunov function induced by a (more common) distance, 
that is Vtcr = Ik — 9*1^2. However, straightforward computations (not presented here) show 
that this candidate Lyapunov function satisfies Vter(p) > for some p. This remark also 
applies to the finite dimensional case. 
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V. Convergence analysis for monotone oscillators 

The result of Theorem [T] has a strong imphcation in the case of monotone PRC's. It imphes 
that the Lyapunov function (I24p has a monotone time evolution if the PRC is monotone. 
In this situation, the Lyapunov function either converges to a lower bound or to an upper 
bound. These two bounds correspond to the two particular behaviors which characterize the 
dichotomy highlighted in Section IIIII They are given by 

< V(p) < + = Q(l) - Q(0) + Q*{1) - Q*{0) = An. 

At the lower bound, the Lyapunov function is equal to zero if and only if the density 
corresponds to the asynchronous (stationary) density (1151) . On the other hand, we will show 
that the function tends to the upper bound 47r if the density p tends to a Dirac function 
(synchronization) . 

A. Exponential convergence to the asynchronous state 

Theorem [T] will be used to study convergence to the asynchronous state for monotone 
decreasing functions K Z{0). However, in order to apply ( l25l) along the solutions, we need 
to show independently that the flux Jo(^) remains strictly positive and uniformly bounded for 
all time. This condition will restrict the set of admissible initial conditions. In particular, the 
initial conditions have to ensure that < Jo < oo when any oscillator crosses 9 = 211 for the 
first time. Formally, we consider the characteristic curves Ag(t) defined by = f(Ae(t),t), 
A,(0) = ^G [0,27r], Ae(te) = 27r. 

The (strictly positive) initial density p(^, 0) > must be such that the value of the flux 
at the intersection of the characteristic curves with = 27r is strictly positive and bounded, 
that is 

p{9, 0) > ^ < j(Ae(te) = 27r,te) = Jo^fe) < oo G [0, 27r] . (44) 

Condition fl44p is the condition that is imposed on the initial conditions to ensure that the 
flux satisfies < Joif) < oo for t G [0,te=o]- But if the flux Jo(t) is bounded for t G [0,te=o]5 
then the flux Jo(^) is uniformly bounded for all time and Theorem [T] can be applied to 
prove the exponential decreasing of the Lyapunov function. The result is summarized in the 
following proposition. 
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Proposition 2. Consider the transport PDF l[T^) - l[T3\) and assume that Z{9) is such that 
(i) the stationary density [W\) exists, (ii) K Z'{9) < 0\/6 E [0, 27r], and (Hi) either Z"[9) > 
e [0,27r] or Z"{9) < W9 e [0,2n]. Then all solutions p{9,t) e C°([0,27r] x M+;M+), 
with WpWli = 1 CL^d with an initial condition satisfying exponentially converge to the 

asynchronous state and the Lyapunov function (f^l) is exponentially decreasing along them. 



Proof: We proceed in two steps. 
Step 1. Let [ii = 0,ti], . . . , [ti,ti], with tj = be tlie successive time intervals so that 
the characteristic curves ^it.j,] satisfy A^ j_^{tj) = and = 27r. We will show that 

Jo{t) is uniformly bounded by 

J"'"' < Mt) < J"'''" Vt > , (45) 

with the bounds J™™ = min^g^ ^o(^) and J™^^ < max^g^ ^o(^)- 

We consider a characteristic curve A^^^(t) with [t, t] C (i.e. t G [tj„i,tj_i] and 

t G [U,ti\). Solving the total derivative equation (HM on A[^j](t) yields 

p{0,t)exp(^- l^'jo{t)KZ'{A^^j^{t))dt^ =p(27rj). (46) 

Next, using ([7j) and expressing the integral in the space variable A^j^ along the characteristic 
curve lead to 

^o(^) f f^^ Jo ^ _ ^o(t) 

uj^KZ{^)m''^^\ k c^ + /^Z(A^^)Jo ^^'^^"^ [^'V "^ + /^^(27r)Jo(t)- 

(47) 

Since p(6', 0) > 0, the exponential evolution (H^ of the density along (all) the characteristic 
curves implies that the density remains strictly positive for alH > 0. Hence, the flux is also 
strictly positive for alH > 0. Then, we define = max^gj^^^ ^o(^) > and we can define 



d = max If ° — —K Z\A^-^) dA^j^ 

\t-t\c[t,_„u] \ Jo u + K Z (Afj ji) Jo ^' ' 

(48) 

' " TV—. K Z'(A,,j.)dAr,T. 

c^ + A'Z(A[,,^) ^ ^'^^ ^'^ 

Since K Z' < \/9, one has Ci = 1 only if Jo{t) is constant on at least one interval [t,t\. In 
this case, the solution has reached the steady state and it follows that 

Jo{t) = J* yte [ti, oo) . (49) 
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Otherwise, one has Cj < 1 and, using (H5i) and computing the integral in equation (H7I) . we 
obtain 



Tmax ^ ^« 



w + /s: Z(27r) Jo(t) - + is: Z(0) Jo{t) yu + K Z{2n) J^^ ' ' ^ ' 

Since Jo(t) < J|^^|]'' by definition, it follows that 

^ 777 — T IIISjX 



u + A' Z(27r) Jo{t) - u + K Z{2ti) J^^ \t^G\t^_,M yto + K Z(0) J^S^ 



<i 

which implies that there exists a constant C™^^ < 1 (see Remark [T]) such that 

Jo(t)<CrV[-f WG[i„t,]. (51) 

It follows from that 

Jo(t) < CrV— VtG[i„t.] (52) 
and a straightforward induction argument using (HUj) and (132]) implies that 

Jo(t) < f fl Cfc"^) ^[of] < -^"^'^ Vt > , 

\fc=2 / 

where some constants C™^^ are equal to one if fl49p is satisfied at time t. The proof for the 
lower bound J™™ in (1451) follows on similar lines. 

Step 2: Since p{0,t) > for all t > 0, the result of Theorem [T] is applied and yields 

J^in min (KZ'{e))v{p)<V{p)<J^,. max (k Z'{e))vip) 

or, equivalently, 

V(p)<0 



^ t^max max 

6»G[0,27r] 



KZ'{9) V(p) < V(p) < - Jmin min Z'(^) 

0e[O,27r; 



since K Z' < 0. The initial condition (1H|) implies that J™™ > and J™^^ < c)0 and the 
Lyapunov function is exponentially decreasing, which concludes the proof. 

Remark 1. When KZ {27v) < 0, in the particular case J^j^ = uj/\KZ{2'k)\, does not 
hold. However, since p(6', t) > for all t > 0, the value uj/\KZ{2'k) \ can be used as the upper 
bound J'^'^^ in (H5]) . 
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Proposition[2]is a strong result showing that, provided that the function K Z{6) is decreas- 
ing, the solution p{-,t) remains, for all time, in a particular set of functions {p|V(p) < C}, 
with the constant C > 0. Inside this set, the solution eventually converges at exponential 
rate toward the stationary solution p*, corresponding to V(p*) = 0. 



The restriction (1441) on the initial condition is rather weak for decreasing functions K Z{6), 
as shown by the following proposition. 



Proposition 3. Assume that K Z'{6) < \/6 e [0,27r]. If K Z{2'k) < 0, then is always 
satisfied. If K Z{2'k) > , then is satisfied if 

P(^,0)<-^^ V^G[0,27r]. (53) 

Proof: Case KZ(27r) < 0. Using the boundary condition (fT3|) . the condition < 
Jo{te) < C)0 is turned into a condition on p{27i,tg) and is equivalent to 

p{9,0)>0 =^ p{2n,tg)>0 y9e[0,2n]. (54) 



For 9 = 27r, condition (IMI) is satisfied, since t2n = 0. Next, we proceed by induction on 9: 
given 9 and assuming that ( l54l) is satisfied for all 9 > 9, we prove that (l54l) also holds at 
9. The total derivative equation (1141) is well-defined on the characteristic curve Ag{t) since 
< Jo(^) < oo for all t <tg. Solving (IT^ along the characteristic curve yields 

p{27T,te) = p{9,0) exp J^' Ut) K Z\Xe{t)) dt^ (55) 

and p{0,0) > implies p{2'K,te) > 0. Condition (Hll) . equivalent to (|M|) . is then always 
satisfied. 

Case KZ(27r) > 0. Using f[T^ . the condition < Jo(^e) < oo is turned into a condition 
on p(27r, te) and (jH]) is equivalent to 

p(^,0)>0 ^ 0<p(27r,t,)<-^^^ V^G[0,27r]. (56) 

The strict condition p(^, 0) > always implies p{27i,te) > 0, as in the case K Z{2tc) < 0. 
Hence, we focus on the additional upper bound on the density p(27r, tg). For 9 = 2n, condition 
implies (13^ . since t2n = 0. Next, we proceed by induction on 9: given 9 and assuming 
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that (jSHD is satisfied for all 6 > 6, we prove that (IHB]) also holds at 9 (provided that condition 
f l55|l is satisfied). Using flSSjl with condition flS^ leads to 

P(2vrj,) < ]^^exp jj' Jo(t) Z'(A,(t)) rft^ . 

Expressing the integral in the space variable along the characteristic curve yields 



Since K Z' is negative, the flux Jq can be replaced by its maximal value, that is, Jq ^ oo 
and the above inequality leads to 

"f^^-*"' < KwrA- 1'' ^)''''^''''^') ■ 

The relation 0571) is well-defined since Z{9) ^ for all 9 G [0,27r]. Finally, computing the 
integral in implies that p{2'7i,tg) < \/[K Z{2ti)]. Condition equivalent to fl5B]) . is 
then satisfied. This concludes the proof. ■ 
When the conditions of Proposition [3] are not satisfied, condition 0441) may fail to hold, in 
which case the flux blows-up (finite escape time to infinity). In the case K Z{2ti) > 0, there 
is a positive feedback between the flux and the velocity: a high value of Jq increases the 
velocity through the coupling, which in turn increases the flux. When the density approaches 
the critical value 

the flux is high enough to blow-up through the positive feedback. 

This finite escape time phenomenon is related to the absorption phenomenon observed for 
finite populations. If two oscillators are close enough, the firing of one oscillator can trigger 
the instantaneous firing of the second: the latter is absorbed by the former. In particular, 
condition flS51) is equivalent for finite populations to an initial condition imposing a minimal 
state distance K/N between any two oscillators, a condition that prevents the absorption 
phenomenon. 

B. Finite time convergence to a synchronous state 

For increasing functions K Z{9), Theorem [T] implies that the Lyapunov function fl2^ is 
strictly increasing, and a synchronous behavior is observed in finite time, for any initial 
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condition. Either the flux Jo(t) becomes infinite in finite time or the density p{0,t) becomes 
infinite in finite time. 

The finite time convergence to synchronization is estabhshed in Proposition HI As a 
prehminary to this result, we need the following lemma. 



Lemma 1. The Lyapunov function satisfies 

V = ||g - g*||Li < 4:71 -2qmin 

with 

g™„ = min min q{ip), min q*{ip)] . 

Proof: The proof of Lemma [1] can be found in Appendix. ■ 
Through Theorem [T] and the preceding lemma, the following proposition establishes the 
finite time convergence to the synchronous state. 

Proposition 4. Consider the transport PDF l[T^) - l[T3\) and assume that Z{9) is such that 
(i) the stationary density [W\l exists, (ii) K Z'{9) > G [0,27r], and (Hi) either Z" {9) > 

G [0,27r] or Z"{e) < G [0,27r]. Then all solutions p{9,t) G C°([0,27r] x [0, 
tfin < oo, with IIpIIli = 1 o,nd with an initial condition p{6, 0) > 0, converge in finite time to a 
synchronous state. That is, if K Z{0) > 0, the flux satisfies J{0,tfin) = oo, or if K Z{0) < 0, 
the density satisfies p(0,t^„) = oo. 

Proof: An infinite fiux Jo{t) or an infinite density p{0,t) is obtained when the density 
p(27r, t) reaches a critical value. If K Z{0) > 0, the value K Z{2'n') is positive since i^Z is 
increasing. Hence, the fiux Jq becomes infinite when the density p(27r, t) exceeds the critical 
value fISS]) . If K Z{0) < 0, the velocity ffTU]) at 6* = is equal to zero when the fiux reaches 
the value Jo(t) = u/\K Z{0) \ or equivalently, given (IT^ . when the density reaches the value 

^P--"= /fZ(2.)-A-Z(0) - 
With a velocity equal to zero at 6' = 0, the relationship ([7]) implies that the density is infinite 
at 6 = 0. (If K Z{0) < along with K Z{27r) > 0, the reader will notice that the value 
has no importance, since (l58l) is greater than (l59l) .) 

Next, we show that the density p(2-7r, t) must necessarily reach the critical value 
or (l59l) in finite time tgn. Let us consider a characteristic curve A(t), with A(t) = and 
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A(t) = 2tt and assume that the synchronous state is not reached within [t,t\, so that p{6, t) G 
C''([0, 27r] X [t, t];]RQ"). (The initial condition p{6,0) > imphes that the density remains 
strictly positive on (all) the characteristic curves.) Applying Theorem [T] and integrating (I25p . 
one has 

min (KZ'ie)) f Mt)dt< l-dV < max (KZ'i9) \ f' Mt) dt . 

ee[o,2n] \ J Jt Jv{t) V ee[o,27r] V / 4 

Since [t, t] is the time interval corresponding to the complete evolution of an oscillator from 
^ = to 6' = 2-7r , the integral of the flux Jq in the above equation is equal to one and it 
follows that 

expf min ( K Z' (6))] < < exp I max (KZ'(e) 

The condition K Z' > implies that the Lyapunov function strictly increases within the 
time interval [t,t\. Considering n successive intervals [0,ti], . . . , [t^jtri], with ti_^_i = ti, one 
obtains 

> exp |n min (kZ'(9) 
Hence, a given value V > V(0) is reached within at most nmax time intervals, with 



log 1 


(V/V(0)) 


1 


mineg[o,27r] 1 


(KZ'i9)) 



Since the time length of each interval [t^j^fc] (A; = 1, . . . , n^ax) is finite, any value V < 47r is 
reached in finite time. By Lemma [H V = V implies that gmin < (47r — V)/2. When considering 
values V close to Air, the minimum of the quantile density q reaches in finite time any given 
value close to zero. Given fl22]) . this implies that the maximum of the density pM reaches 
in finite time any given value (provided that the critical value ( 1351) or (15^ is not already 
reached). In particular, the value 

pM = Pcexp ^^max^ (^K Z'{e)j^ , (60) 

with pc denoting the critical value (I58p or (I59j) . is reached in finite time. Then, the maximum 
value pm (obtained at 9m at time t^) decreases along the characteristic curve A(t), with 
A(tj\/) = 6m- Solving the total derivative equation (IT^ . one shows that the variation along 
the characteristic curve is bounded: 

p(27r, tfi,) = PM exp y^*'" Jo{t) K Z'[A{t))d?j > pm exp max^ [k Z'(e))^ , 
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where the inequahty is obtained since the integral of Jq on tgn] is less than one. With 
pM given by fl^ . the density exceeds the critical value in finite time tgn, which concludes 
the proof. ■ 
Proposition H] shows that oscillators with monotone increasing functions K Z converge to 
a synchronous behavior corresponding to a solution with a finite escape time to infinity. If 
the coupling is excitatory (0 < K Z{0) < K Z{27r)), the fiux Jq blows-up and is infinite at 
time tfin, as well as the velocity t'(27r, tgn). If the coupling is inhibitory [K Z{0) < 0), the 
velocity v{0, t) reaches the value zero at time tgn and the oscillators accumulate at 6* = 0, 
yielding an infinite density p(0,tfin)- 

VI. Implications for populations of oscillators 

The results of Section IIVI and Section |V] are technical. In this section, their implications 
on the behavior of (monotone) pulse-coupled oscillators are discussed with more details. 

Emphasis is put on the strong parallel between the analysis of the present paper on infinite 
populations and earlier results on finite populations. We also consider more general models 
and discuss the relevance of an instantaneous impulsive coupling. 

A. Parallel with finite populations 

The results of Proposition [2] and Proposition H] apply to oscillators characterized by a 
monotone PRC. In the case of integrate-and-fire oscillators x = F{x), the results thereby 
apply to oscillators with a monotone vector field F (including the popular LIF oscillators) 
and reinforce the numerical observations presented in Section IIII-AI 

When K Z' < [KdF/dx > 0), the oscillators converge to the asynchronous state. 

Theorem 2. Consider a continuum of identical pulse-coupled integrate-and-fire oscillators 
X = F{x) characterized by (i) a monotone dynamics KdF/dx > and (ii) a PRC with a 
curvature of constant sign. Then, provided that the asynchronous state exits (cf. Corollary 
[ip, either the initial conditions do not fulfill (IHI) and the flux tends to infinity in finite time, 
or the continuum exponentially converges to the asynchronous state. 

When K Z' > {) [KdF/dx < 0), the oscillators achieve synchrony. 



January 12, 2013 



DRAFT 



28 



Theorem 3. Consider a continuum of identical pulse-coupled integrate- and- fire oscillators 
X = F{x) characterized by (i) a monotone dynamics KdF/dx < and (ii) a PRC with a 
curvature of constant sign. Then, provided that the asynchronous state exits (cf. Corollary 
[IP, the continuum converges to the synchronous state in finite time. 

The curvature condition on the PRC is verified for the LIF oscillators (with dynamics 
X = S — x). Indeed, the PRC is given by 



and satisfies Z" > 0. Hence, Theorem [2] and Theorem [3] prove the global dichotomic behavior 
of the continuum of pulse-coupled LIF oscillators: the oscillators converge to the synchronous 
state when the coupling is excitatory {K > 0) and to the asynchronous state when the 
coupling is inhibitory {K < 0). These global results complement the local results on the 
continuum model presented in P, [T7|, [53] . 

The results of Theorem [2] and Theorem [3], for infinite populations, are the exact analogs of 
the results presented in our previous study [12] and in [20], for finite populations. They prove 
that the behaviors of monotone oscillators are dichotomic, not only for finite populations 
but also for infinite populations, thereby drawing a strong parallel between the analysis of 
finite populations and the analysis of infinite populations. 

B. Application to other oscillators 

Beyond the case of integrate-and-fire oscillators, the results developed in the present paper 
apply to more general phase dynamics. The PRC can be (numerically) computed from 
general high-dimensional state models possessing a stable limit cycle, in which case the phase 
dynamics (jlj) represents a one- dimensional reduced model valid in the neighborhood of the 
limit cycle [30], [M]- However, two limitations appear, which require (i) a weak coupling 
strength and (ii) a monotone PRC. 

(i) Phase reduction of multidimensional models (other than integrate-and-fire models) is 
valid as long as the coupling is weak. Since the derivation of the transport equation relies 
on the phase dynamics, it is relevant only for a coupling strength K <^1. 

(ii) The results only apply to strictly monotone PRC's. As a consequence of the mono- 
tonicity, the PRC must be characterized by Z{0) ^ Z{27t) and cannot be periodic. However, 




with 
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a PRC computed on a limit cycle is periodic, unless the limit cycle is discontinuous. For 
instance, such a discontinuity is artificially created in the one-dimensional integrate-and- 
fire model (or in the Izhikhevich model [H], [13]). A discontinuous limit cycle can also be 
obtained for models with separated time scales (e.g. relaxation models such as FitzHugh- 
Nagumo oscillators [H], [22] or Van der Pol oscillators, spiking oscillators). Since the time 
spent on a part of the limit cycle is negligible with respect to the time spent on the rest 
of the cycle, the fast part can be replaced in good approximation by a discontinuity. The 
oscillators are thereby characterized in good approximation by a non periodic PRC, that 
can be monotone. 

This situation is illustrated in the following example, for oscillators characterized by a 
limit cycle close to a homoclinic bifurcation. 

Example. In [5], limit cycles close to a homoclinic bifurcation are approximated by a dis- 
continuous limit cycle, to which corresponds a monotone (discontinuous) PRC that satisfies 
the additional curvature condition, so that the results apply. 

A homoclinic bifurcation occurs when there exists, for a given parameter value, a homo- 
clinic orbit to a saddle point with real eigenvalues. At the bifurcation, a limit cycle appears. 
Assuming that there is a single unstable eigenvalue such that < |AsJ, with Ag^. being 
the stable eigenvalues, the limit cycle is stable [TT]. Since the trajectory is much slower near 
the saddle point, the limit cycle is discontinuous in good approximation and the PRC is 
approximatively given by the discontinuous function (see [5] for more details) 



where C > is a mo del- dependent constant. 

The PRC f l6ip is monotone decreasing and has a positive curvature. It follows that the 
oscillators close to a homoclinic bifurcation and interacting through a weak excitatory 
impulsive coupling (0 < A' ^ 1) satisfy the hypothesis of Proposition [21 they exponen- 
tially converge toward the asynchronous state, provided that the asynchronous state exists. 
The condition flTHl) of Proposition [T] is always satisfied with a weak coupling, so that the 
asynchronous state always exists. For a weak inhibitory coupling (— 1 ^ < 0), Proposition 
m implies that the oscillators reach a synchronous state in finite time. 




(61) 



January 12, 2013 



DRAFT 



30 



In the popular Morris-Lecar model [21] , which is among the most widely used conductance- 
based models in computational neuroscience, a homoclinic bifurcation can occur for low 
external currents [25], [21] • The above results apply in this situation. o 

VII. Conclusion 

In this paper, we have studied the global behavior of infinite populations of monotone 
pulse-coupled oscillators. The behavior of the oscillators is dichotomic: either the oscilla- 
tors achieve perfect synchrony (synchronous state) or the oscillators uniformly spread over 
S'"'^(0,27r) (asynchronous state). 

The infinite population is represented by a continuous density. In this framework, a 
necessary and sufficient condition ensures existence and uniqueness of the stationary density 
(asynchronous state). 

For the global convergence analysis of the (nonlinear) transport equation of the density, 
we propose a Lyapunov function that is the total variation distance between quantile density 
functions. The Lyapunov function has two extreme values that correspond to the two steady- 
state behaviors of the system (the synchronous state and the asynchronous state). The 
stability results obtained for general phase oscillators are applied to particular models of 
importance (e.g. leaky integrate-and-fire model). 

The main result of the paper stresses the importance of a distance (the total variation 
distance) to analyze a transport PDE under monotonicity assumptions on the PRC. How- 
ever, the time evolution of the proposed Lyapunov function is no longer monotone when the 
monotonicity assumptions fail, even though the observed dichotomic behavior seems more 
general. This restriction raises interesting open questions about the generalization of the 
Lyapunov function and about the use of distances for a larger class of transport PDE's. 

Appendix 

The Lyapunov function is written as 

V= / \q-q*\dip =f {q-q*)dip+ f (q* - q) dip 
Jo Ja+ Ja- 

< {Q- Qmin) d(p+ {q* - gmin) dif , (62) 

ja+ ja- 
with ^+ = e [0, l]|g - g* > 0} and A" = e [0, l]|g - q* < 0}. 
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Moreover, one has 



A+ 



qd(p + 




qmmd(p< [ qd(p = Q{1) - Q{0) = 271 , 



(63) 



and 



JA~ Ja+ Jo 



[ q*dip+ [ q^,i,,dip< [ q* dip = Q*{1)-Q* {0) = 271 



(64) 



Next, injecting (1^ and (IMI) in inequality yields 

V < 27r — 2 / gmin dip + 27r — 2 / gmin (iy? = 47r — 2 g, 
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